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We propose a method of measuring diffusion of the rela- 
tive phase between two Bose-Einstein condensates occupying 
different nuclear or spin hyperfine states coupled by a two- 
photon transition via an intermediate level. Due to the macro- 
scopic quantum coherence the condensates can be decoupled 
from the electromagnetic fields. The rate of decoherence and 
the phase collapse may be determined from the occupation of 
the intermediate level or the absorption of radiation. 
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In this paper we propose a method of measuring the 
diffusion of the relative phase between two Bose-Einstein 
condensates (BECs) formed in different hyperfine states 
of the same atom. These two states are coupled by a two- 
photon transition via an intermediate atomic level. Due 
to the macroscopic quantum coherence of the BECs the 
two-photon transition between the two hyperfine states is 
suppressed by quantum interference effects. This effect is 
similar in origin to that occuring in electromagnetically 
induced transparency (EIT) and lasing without in- 
version Atomic interactions give rise to phase diffu- 
sion which destroy the coherence, so that the two-photon 
transition is no longer completely suppressed. The rate 
of phase diffusion may be determined by monitoring the 
population in the intermediate level or the absorption 
rate. 

Since the first realizations of dilute gas alkali BECs ^] 
the experiments have been broadened to include two and 
multiple condensate systems [||-||. Myatt et al. Q pro- 
duced two overlapping BECs of |F=1, m=-l) and |F=2, 
m= 2) states of *^Rb using sympathetic cooling. The sta- 
bility of this pair is due to an unexpectedly small inelastic 
collision rate between these states . Recently, another 
BEC pair of |1, -1) and |2, 1) states of ^''Rb has been re- 
alized at JILA 10-^. These two states have essentially 
identical magnetic moments and an adjustable spatial 
overlap allowing the creation of a fully interpenetrating 
binary mixture. An especially interesting property is that 
these two BECs can be coupled by a two-photon transi- 
tion (one microwave and one radiofrequency photon) . 

In general, inelastic collisions will limit the possibili- 
ties of magnetically trapping BEC pairs. However, opti- 
cal dipole traps, which use optical forces to trap atoms, 
have a major advantage over magnetic traps since they 
can stably trap atoms in arbitrary hyperfine states. An 
evaporatively cooled ^"^Na gas has been succesfuUy con- 
fined in a dipole trap with a simultaneous observation of 



BECs in several different hyperfine states [|j . 

A fascinating property of BECs is that they exhibit a 
macroscopic quantum coherence that is absent in ther- 
mal atomic ensembles [0. Since one needs a phase ref- 
erence to observe a phase, binary mixtures of BECs are 
especially useful in the studies of coherence properties. 
The atom-atom interactions in finite-sized BECs affect 
the matter wave coherence. The width of the number 
distribution in the ground state has a dispersive effect on 
the BEC self-interactions and the relative phase under- 
goes quantum collapses and revivals Additional 
sources of phase diffusion are spatial mode fluctuations 
[Q and finite temperature decoherence due to the in- 
teractions between condensate and noncondensate atoms 

In this paper we consider a system closely related to 
the recent experiments of overlapping BECs in different 
hyperfine levels coupled by a two-photon transition ■ 
Analogous BEC pairs could possibly be produced also in 
dipole traps 1^. As a consequence of the macroscopic 
quantum coherence of BECs the two hyperfine levels cou- 
pled by a two-photon transition exhibit two-photon co- 
herence. By adjusting the initial conditions of the BECs 
and the driving electromagnetic (EM) fields the atoms in 
the BECs can be decoupled from the EM fields. However, 
as a consequence of the decoherence and the collapse of 
the relative phase the two-photon coherence is reduced, 
and the cw EM fields start inducing atomic transitions. 
This provides an excellent scheme for measuring phase 
dynamics in the present experimental set-ups . The 
two-photon coherence of two optically coupled BECs has 
been previously predicted to result in various dramatic 
properties of the scattered light ||l^,|l8| . Measurements of 
magnetic coherence-related phenomena in atomic BECs 
have been previously addressed, e.g., in spin-polarized 
hydrogen ||l9|l and in the electronic spin resonance of al- 
kali gases po| . 

We consider a three-level system with the energies of 
levels |1), 1 2), and |3) denoted by wi, W2, and L03. The mi- 
crowave or rf cw EM fields and induce magnetic 
dipole transitions between the levels 1 <-^- 2 and 2 <-> 3, 
respectively. For simplicity, we assume that the radiative 
lifetimes are long, so that the level widths can be ignored. 
The EM interaction introduces the following terms into 
the Hamiltonian density: 
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'fJ'12 ■ Ba'0IV'2 - M23 ■ BbV^V-s + H.c. , (1) 



where "01 (r) is the field operator for hyperfine state |i) 
and /Xjj denotes the magnetic moment for the transition 
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We consider a situation where two macroscopic, per- 
fectly overlapping BECs occupy levels |1) and |3) and 
level 1 2) is initially empty. In Refs. ||-§ the double BEC 
system in levels |1) and |3) is prepared from the single 
BEC in state |1) by a two-photon transition. If the oc- 
cupation of level 1 2) is small, the collisional interactions 
are mainly between atoms in levels |1) and |3). We write 
the interaction Hamiltonian density in terms of the field 
operators V'i(r) and V'3(r) 

Hint = X! y'^l^J^^'^i + wiai/'iV'IV'ai/'i , (2) 

i=l,3 

where Ui — Anh^ai/m. We could directly obtain the 
equations of motion for the condensate mean fields from 
Eq. (^) , but the interactions are significantly simplified in 
the case of approximately equal scattering lenghts result- 
ing in wi ~ M2 ~ ui2 = u. In ^''Rb the scattering 
lengths satisfy Oi : ai2 : 02 :: 1.03 : 1 : 0.97, where ai (02) 
denotes the intraspecies scattering length for state |1, - 
1) (|2, 1)) and 012 is the interspecies scattering length. 
We approximate the field operators il^iiv) ~ (f>i{r)ai for 
levels |1) and |3) in terms of the BEC annihilation op- 
erator ai and the corresponding spatial wave function 
4>i{r). It is also assumed that 0i = ^3 = 0. Our aim 
is to write the equations of motion for the expectation 
values Uij = {a\a-j)/N, where the total initial number 
of BEC atoms is denoted by N . The effect of noncon- 
densate atoms is treated by a phenomenological damping 
in the equation of motion for the coherence between the 
two BECs. By approximating the scattering lengths to 
be equal and by assuming {a\a2) N we may then ap- 
proximate the interaction Hamiltonian [Eq. (S)] by 



{N^ - N -2Na\a2) 



(3) 



where N — a\ai is assumed to be a constant of 
the motion and Uk — u j (P r\(l)\'^ . We see that for 
perfectly overlapping BECs with approximately equal 
scattering lengths the effect of BEC self-interactions is 
strongly suppressed. Finally, the detunings are defined 
by 821 = oja- ^21 and (S32 = lob - W32, where (wb) 
is the frequency of the field (Bb) and the transition 
frequency between the hyperfine levels 1 <-^- 2 (2 <-> 3) is 
^21 (^32)- In the rotating- wave approximation we then 
obtain the following equations of motion for the expec- 
tation values (T,- 
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dii = fi^Im(cr2i), 

0-22 = -f^Alm(cr2l) + f^slm(cr32), 
0-33 = -^^Blm((T32), 



(4a) 
(4b) 
(4c) 
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-i((52i -f- NK)a2i 



2 ((T22-CT11) ^^31, (4d) 



^^32 = - Ahi - NK)a-i2 + ^-(^33 - CT22) + -^^31, (4e) 

J. . , i^A ,,r\ 

0-31 = -i((>32 + O21 - ll)<J3l + ^~^32 2"^^^^' ^ ' 



where the Rabi frequencies are given by Qa = 
2 Jrf^r(/)2^i/i2i-B^/?i, Ob = 2 j (Pr(j)2(l)3li2z-'^B I'h, and 
Im denotes the imaginary part. For simplicity, in Eq. (Q) 
we have set Q,a and JIb to be real. Here the equations 
correspond to the cascade or ladder three-level system 
with cji < [x)2 < ^3- In the case of A (V) three-level 
scheme the sign of ^32 {52i) should be changed. The dom- 
inant mean-field contribution of the BEC self-interactions 
is to shift the resonance conditions in Eqs. (^) and (^). 

Level 1 2) is assumed to be initially empty and two 
BECs occupy levels |1) and |3). The off-diagonal element 
CT31 describes the macroscopic coherence between the two 
BECs. As explained earlier this collapses and decoheres 
due to the atom-atom interactions. We have included the 
effect of the decay of the matter wave coherence in Eq. (^) 
in terms of a phenomenological damping parameter 7 in 
the equation of motion for CT31. This damping parame- 
ter includes contributions from both the quantum effects 
of BEC self-interactions and collisions between conden- 
sate and noncondensate atoms. We assume that these 
dominate over other damping mechanisms as long as the 
population in level |2) remains small. 

For an initial condition for Eq. (^ we set cth = 1733 = 
CTsi = 1/2 indicating a well-established coherence be- 
tween the two BECs with a vanishing relative phase. 
This corresponds, e.g., to a situation where a BEC is 
first prepared in level |1) and half of the BEC atoms 
are then coherently transferred to level |3), so that the 
atoms remain entangled. Initially there are no atoms in 
the intermediate level (T22 = 0. The EM fields are two- 
photon resonant {821 — —5^2) and U,a — — f^B = It 
is easy to see from Eq. (||) that for 7 = this corre- 
sponds to a steady-state situation. Both BECs are de- 
coupled from the EM fields and the absorption described 
by Im(cr2i) and Im(CT32) vanishes |^^. However, for non- 
zero 7, Im((T2i) and Im((732) also become non-zero, and 
atoms start accumulating in level |2) due to the absorpion 
of EM radiation. The phase damping parameter 7 may 
be determined by measuring the oscillations of the EM 
fields or the population in level |2). 

Equations (^) may be integrated numerically, but it 
is also illuminating to look at analytic estimates. We 
consider a situation where the EM fields are resonant, 
821 + Nk — 632 — Nk = 0, and the damping 7 is much 
smaller than the Rabi frequency fl. We look for an expo- 
nential solution to Eq. (^) to leading order in the small 
parameter ^/ft. By also determining the coefficients to 
first order in j/il we obtain (731 (i), describing the matter 
wave coherence between the BECs, as 



^731(0 



1 e-37t/4 , 



The absorption of the field B^ is proportional to 

Im[a2i(t)] ^ + ^ e-^*/« cos (^2^) . (6) 

The absorption results in oscillating EM fields with the 
amplitude of the oscillations given by7/(8f])e-''*/8. The 
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real parts satisfy Re[cr2i(t)] 
occupation in level |2) is 



=Re[cr32(t)] = 0. Finally, the 
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-7*/8 



sin {V2nt). (7) 



Due to the decoherence of the BECs the EM fields absorb 
radiation and atoms start occupying level |2). By mea- 
suring the number of atoms in state |2) at time t after 
switching on the driving EM fields, one could determine 
the damping rate of the matter wave coherence 7. Al- 
ternatively, the damping rate could be observed from the 
amplitude of the oscillating EM signal. 

In Fig. |l| we have plotted one example of the signal 
corresponding to a particular value = 2O7. We plot 
Im[cr2i(^)l and (722 (0 obtained by numerically integrat- 
ing Eq. (Q) for 621 + Nk = 632 - Nk = 0, and for the 
initial condition (Tii(0) = cr33(0) — cr3i(0) — 1/2 and 
o'2i(0) = (732 (0) — 0. The oscillating signal (a) and the 
accumulating population in the intermediate level (b) are 
clearly observed. The graphs are also well represented by 
the approximate analytic results, Eqs. (^) and (|^). 

vin-^ Fig.la 
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FIG. 1. A particular example response of the medium in 
the case of f2 = 2O7. (a) Im[(T2i(t)] proportional to the ab- 
sorption of the EM radiation from the field driving the tran- 
sition between levels |1) and |2). (b) 0-22 (i) proportional to 
the population of the intermediate level. The EM fields are 
assumed to be resonant. 

If the fields are off-resonant from intermediate level 
1 2) the absorption of the radiation and the occupation of 



state |2) are reduced, but Re((T2i) and Re((T32) are in this 
case non-zero. 

In EIT the interference between the off-diagonal den- 
sity matrix elements leads to an initially opaque medium 
being rendered almost transparent |^,^. The transition 
amplitudes driven by two oppositely phased EM fields 
destructively interfere. The effect of the medium on the 
EM fields is canceled. The crucial quantities for the in- 
terference are the coherences (T21 and (732 describing the 
transition dipole matrix elements excited by the driving 
EM fields. In EIT the coherence between levels |1) and 
1 3) is present only as a consequence of the EM coupling 
via level |2). For BECs in hyperfine states |1) and |3) 
the coherence 1731 is present from the start without be- 
ing created by the driving fields due to the macroscopic 
matter wave coherence. As a result, the EM fields couple 
(T31 to the oscillating dipoles. The presence of the macro- 
scopic quantum coherence of BECs can then completely 
inhibit the EM fields from establishing any coherence be- 
tween levels |1) and |2) or |2) and |3). On the other hand, 
the rate of which coherences tT2i and (732 are induced de- 
scribes the decoherence rate of BECs. 

In recent experiments Hall et al. B studied phase dif- 
fusion in a binary BEC mixture of Rb. A BEC was 
first prepared in level |1, -1) and a part of the conden- 
sate was then transferred to level |2, 1). The relative 
phase between the two separated halves was determined 
by interfering the atoms at a later time. The phase diffu- 
sion rate was estimated by varying the evolution time of 
the two BECs before the interference measurement. In 
every interference measurement the BECs were destruc- 
tively imaged and the repetitions of independent runs 
produced information about the uncertainty of the phase. 
Only weak phase diffusion was observed. 

In addition to the enviromentally-induced decoherence 
due to the interactions between condensate and thermal 
noncondensate atoms ||l5| , [l6|] the quantum collapse due to 
the BEC self-interactions has an important effect on the 
phase diffusion |jl^,|l3|. This rate dramatically depends 
on the relative strength of the three scattering lengths 
ai, as, and 013 in Eq. (H). Under conditions where the 
scattering lengths are equal, the two BECs are perfectly 
overlapping, and only levels |1) and |3) are occupied, the 
interaction Hamiltonian in Eq. depends only on the 
constant total atom number. This suppresses the phase 
collapse. 

Recent experiments have realized overlapping BECs in 
different hyperfine states j^-^. A two-photon transition 
between a BEC pair of |1, -1) and |2 1) states of ^'''Rb 
via level |2, 0) has been implemented The interme- 

diate state |2, 0) of ^"^Rb for the two-photon coupling is 
untrapable. The atoms in |2, 0) can escape the trap. This 
would correspond in our scheme to an additional damp- 
ing F for level |2). If the population in level |2) is small, 
this damping would be approximately independent of the 
number of atoms in |2). In that case we could add the 
following additional term to Eq. (4b): cr22 — . . . — Ta22- 
Then the phase diffusion could possibly be measured by 
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monitoring the number of atoms escaped through level 
1 2, 0) or by counting the atoms remained in levels |1, - 
1) and |2, 1). Optical dipole traps can stably trap 
atoms in arbitrary hyperfine levels. Suitable BEG pairs, 
without losses of atoms, could possibly be produced in 
dipole traps to implement the proposed scheme for the 
measurement of the phase diffusion. 

In conclusion, we have proposed a method of measur- 
ing the "phase memory" of a BEG pair. This method 
relies on the quantum interference of transition ampli- 
tudes and is similar in origin to that occuring in EIT and 
lasing without inversion. Unlike the previous measure- 
ments of phase diffusion |^ , our model allows continuous 
and nondestructive monitoring of the phase dynamics. 
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